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TTTOLO 

Metodo e dispositivo per la stima della funzione di trasferimento di un canale di 
trasmissione 

DBSCRIZIONE 

5 La presente invenzione si riferisce ai sistemi di telecornunicazioni e in 

particolare riguarda un metodo e un dispositivo per la stima della funzione di 
trasferimento di un canale di trasmissione. 

Come b noto, nei sistemi di telecornunicazioni radiomobili una delle tecniche pid 
utilizzate b la tecnica di accesso W-CDMA (Wideband - Code Division Multiple 
10 Access) in virtii della sua elevata efficienza spettrale nei confronti di altre tecniche di 
accesso multiplo. Neirambito di questa tecnica, e in particolare nella modality FDD 
(Frequency Division Duplexing), sono state adottate diverse soluzioni in grado di 
aumentare la capacita del sistema, fra cui b da citare la rivelazione coerente del segnale 
ricevuto dalla stazione radio base o dal terminale mobile. 

15 Qsssto ripo di rivelazione richiede un'accuratissima stima della funzione di 

trasfer; »*. canale di trasmissione, denominata nei seguito "stima del canale", in 
present o>anescenze e riflessioni multiple, dovute ai movimento deH'utente 
nell'area di copertura, senza peraltro avere informazioni della sua velocity di 
spostamento. 

20 La stima del canale, come altre funzioni del sistema, possono essere assegnate 

ad elaboratori (per esempio DSP, microcontrollori, ecc.)» nei qual caso esse vengono 
realizzate in modalita software, oppure possono essere assegnate a dispositivi dedicati 
(per esempio FPGA, ASIC, ecc), e quindi la modality b prevalentemente hardware. 
Compito del progettista, nei corso della definizione deH'architettura del sistema di 

25 telecornunicazioni, e la ripartizione ottimale dei compiti fra i dispositivi programmati e i 
dispositivi cablati. 

Le diverse architetture realizzabili devono essere analizzate sotto Taspetto del 
carico computazionale delPeventuale elaboratore DSP, del volume di traffico sul bus 
dedicato al trasferimento dei dati, della ripartizione ottimale tra le risorse hardware e 
30 software. 
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Considerando per esempio un sistema complesso come un modem per stazione 
base o per terminate mobile UMTS, ne risulta che il bus di comunicazione tra il DSP e 
le risorse hardware rappresenta il collo di bottiglia del sistema. Ad esempio, nel caso di 
stazione radio base, il carico del bus si fa sempre piil grande con Paumento del numero 
5 degli utenti e con la presenza di funzionalitk che richiedano il trasferimento di notevoli 
quantita di dati tra DSP e risorse hardware. 

La stima del canale & una di queste funzionalite, in quanto essa deve essere 
ripetuta per ognuno dei finger (sottocanali) di ciascun utente. 

Consideriamo il caso peggiore in cui stia comunicando in una cella il numero 
10 massimo di utenti Nu = 128 e ciascun utente abbia un numero di finger Nf = 8 per 
ciascun Rake receiver (ricevitore a rastrello, tipico ad esempio delle stazioni base 
UMTS). La quantity di dati da trasferire per ciascun finger in un time slot (intervallo 
temporale), di durata Ts = 666 us, pud essere stimata pari a un numero di bit Nb = 400. 
Da cio deriva un notevole flusso di dati sul bus, pari a: 
15 NbNf Nu/Ts=400-8- 128/666-10-6 = 615 Mbit/s 

Al giorno d'oggi, sono conosciuti diversi algoritmi, in grado di effettuare la 
stima del canale, fra questi, vi sono ad esempio quelli descritti nel documento 
EP 0 912 019, in cui si utilizzano due differenti metodi di interpolazione, lineare e 
mediante filtri di Kalman, quello descritto nel documento EP 1 032 168, in cui 
20 Tinterpolazione viene fatta con il metodo dei polinomi di Lagrange, e quello descritto 
nel documento US 5,886,9 1 1 in cui viene utilizzato il metodo di bisezione. 

Le soluzioni descritte nei suddetti documenti presentano pero alcuni 
inconvenienti: o non sono applicabili direttamente in sistemi UMTS, come il melodo di 
bisezione, o sono eccessivamente complessi per una realizzazione di tipo hardware. 

25 Considerato infatti il notevole flusso di dati gi& presente sul bus di una stazione 

base, sarebbe preferibile non sovraccaricare ulteriormente tale bus, cosa che accadrebbe 
effettuando complesse operazioni di interpolazione neirambito della stima del canale, 
quali ad esempio le operazioni necessarie per effettuare un'interpolazione con il metodo 
dei polinomi di Lagrange. 

30 II dispositivo pub essere realizzato completamente in hardware, e quindi 

facilmente integrato in un Rake receiver per stazione base, mantenendo le stesse 
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prestazioni delle soluzioni software che si basano su processori DSP. Non essendo 
richiesto alcun trasferimento di dati con il processore DSP, il bus di comunicazione 
viene quindi lasciato libero per altri compiti. 

Sono particolare oggetto della presente invenzione un metodo e un dispositivo 
5 per la stima della funzione di trasferimento di un canale di trasmissione, come descritto 
nelle unite rivendicazioni. 

H melodo e il dispositivo, oggetto della presente invenzione, fanno uso di un 
algoritmo di bassa complessita per la stima della funzione di trasferimento di un canale 
di trasmissione, adatto sia per la tratta di trasmissione verso una stazione base, tratta 
10 denominata "up-link", sia per la tratta di trasmissione verso un terminale mobile, tratta 
denominata "down-link". 

Queste ed altre caratteristiche della presente invenzione risulteranno evidenti 
dalla seguente descrizione di una forma preferita di realizzazione della stessa, data a 
titolo di esempio non limitativo, e dai disegni annessi in cui: 
15 la figura 1 e una rappresentazione schematica dei simboli . r!S canale dati^ 

DPDCH e di un canale di controllo DPCCH, nel caso di una tratte \i > ■ ^ - ifesione verso - 
una stazione radio base; 

la figura 2 e un grafico che illustrauna plurality di coefficienti di canale ottenuti 
per interpolazione lineare fra due coefficienti di canale noti; 
20 la figura 3 £ un diagramma di flusso che illustra un primo algoritmo di 

interpolazione realizzato secondo la presente invenzione; 

la figura 4 fe una rappresentazione schematica che illustra una sequenza di 
operazioni per il calcolo di coefficienti di canale mediante 1'algoritmo di figura 3, nel 
caso Npilot = 3; 

25 la figura 5 & una rappresentazione schematica che illustra una sequenza di 

operazioni per il calcolo di coefficienti di canale mediante 1'algoritmo di figura 3, nel 
caso Npilot = 4; 

la figura 6 6 una rappresentazione schematica che illustra una sequenza di 
operazioni per il calcolo di coefficienti di canale mediante 1'algoritmo di figura 3, nel 
30 caso Npilot = 5; 
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la figura 7 e una rappresentazione schematica che illustra una sequenza di 
operazioni per il calcolo di coefficienti di canale mediante ralgoritmo di figura 3, nel 
caso Npilot = 6; 

la figura 8 & una rappresentazione schematica che illustra una sequenza di 
5 operazioni per il calcolo di coefficienti di canale mediante ralgoritmo di figura 3, nel 
caso Npilot = 7; 

la figura 9 e una rappresentazione schematica che illustra- una sequenza di 
operazioni per il calcolo di coefficienti di canale mediante ralgoritmo di figura 3, nel 
caso Npilot = 8; 

10 la figura lOeun diagramma a stati che riassume le operazioni necessarie per il 

calcolo dei coefficienti di canale nei casi illustrati nelle figure da 4 a 9; 

la figura lie una rappresentazione schematica che illustra una sequenza di 
operazioni per il calcolo di coefficienti di canale mediante un secondo algoritmo 
realizzato in accordo con la presente invenzione, nel caso Npilot = 6; e 

15 la figura 12 e uno schema a blocchi k.v circuito che implementa un algoritmo 

di interpolazione realizzato in accordo cfeyj;- ^vu^ate invenzione. 

II metodo ed il dispositivo secondo la presente invenzione verranno ora descritti 
in dettaglio facendo riferimento in particolare al sistema di telecomunicazioni 
radiomobili UMTS (Universal Mobile Telecommunications System), relativamente all 
20 atratta up-link nella modality FDD. 

La rivelazione coerente per Tinterfaccia radio dei sistemi UMTS pud essere 
effettuata utilizzando i simboli pilota che sono multiplati nel tempo nel canale di 
controllo DPCCH (Dedicated Physical Control Channel). 

In figura 1 e illustrata schematicamente la relazione temporale fra i simboli del 
25 canale dati DPDCH (Dedicated Physical Data Channel) ed i simboli del canale di 
controllo DPCCH in un intervallo temporale detto "time slot". Nel caso illustrate in 
figura sono presenti, all'inizio dello slot, tre simboli pilota Xq(0), Xq(1) e Xq(2), seguiti 
da altri simboli utilizzati per il controllo, in particolare quattro simboli TFCI 
(Transport - Format Combination Indicator), due simboli FBI (Feedback Information) 
30 ed un simbolo TPC (Transmit Power Control). 
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I simboli pilota Xq(k), con k=0,l,2,.., Npilot-j (dove Npilot & il numero di 
simboli pilota nello slot) sono noti anche al ricevitore che utilizza tale inforraazione 
unilamente ai simboli ricevuti Z(k) al fine di valutare l'effetto del canale su ciascun 
pilota X Q (k). 

5 Eseguendo il prodotto complesso del simbolo ricevuto Z(k) per il corrispondente 

simbolo pilota Xq(1c) si ottengono i coefficienti di canale; in particolare, le componenti 
in fase ed in quadratura del prodotto corrispondono alle componenti Cq e Ci dei 
coefficienti di canale C(k)=Q(k)+jCQ(k). 

Risulta pertanto evidente che & possibile stimare i coefficienti di canale 
10 solamente in corrispondenza dei simboli pilota trasmessi sul canale DPCCH. Tuttavia, 
poiche i simboli pilota non vengono trasmessi in continuazione ma solamente in 
corrispondenza della parte iniziale di un time slot, per effettuare la stima del canale 
necessario calcolare, mediante una qualche forma di interpolazione, i coefficienti di 
canale corrispondenti ai simboli pilota mancanti. 
15 H numero di simboli pilota all'interno di un singolo slot puo variare da un 

minimo di fr& r v.** Massimo di otto, pertanto il numero di coefficienti di canale che A 
dovranno es.-io.-. LVnati mediante interpolazione sari variabile fra sette e due, essendo" 
dieci il numero totale di simboli presenti in uno slot di un canale di controllo DPCCH. 

In figura 2 e illustrato ad esempio il caso in cui sette coefficienti di canale (in 
20 figura sono rappresentate solamente le componenti in fase) indicati complessivamente 
con il riferirnento 2, sono stati calcolati mediante interpolazione fra Tultimo coefficienle 
Ci(Npilot-1) dello slot corrente, corrispondente ad un valore di ascissa A, ed un primo 
coefficiente dello slot successivo, indicato in figura come Ci(10) e corrispondente ad un 
valore di ascissa B. 1 simboli pilota, e pertanto anche i coefficienti di canale, 
25 corrispondono sempre, per convenzione, a valori di ascissa interi (0, 1, 2, 3, ...) su un 
asse orizzontale temporale (Time). 

II metodo ed il dispositivo di interpolazione realizzato secondo la presente 
invenzione consentono di calcolare le componenti Ci(k) e C Q (k) dei coefficienti di 
canale senza fare uso di operazioni complesse quali moltiplicazioni e divisioni, ma 
30 solamente utilizzando somme e divisioni per due (quest'ultime facilmente 
implementabili in hardware mediante un'operazione di shift a destra di un registro), 
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consentendo di ridurre notevolmente la complessite dell'architettura hardware dell'unil£ 
di interpolazione. 

Verri ora descritto in dettaglio il metodo di interpolazione secondo la presente 
invenzione, con particolare riferimento al diagramma di flusso illustrato in figura 3. 
5 L'algoritmo illustrato in figura 3 consente di calcolare per interpolazione una plurality di 
punti compresi fra un ultimo coefficiente di canale, di ascissa A, di uno slot corrente L, 
ed un primo coefficiente di canale, di ascissa B, di uno slot successivo L+l a detto slot 
corrente. Essendo pari a tre il numero minimo di coefficienti di canale presenti in 
ciascuno slot, il coefficiente, di ascissa A-l, immediatamente precedente all'ultimo 
1 0 coefficiente di canale di ciascuno slot & sempre noto, e pub quindi essere utilizzato per il 
calcolo di interpolazione. 

Airinterno deiralgoritmo di figura 3 vengono utilizzate alcune variabili, X, Y, Z, 
W che hanno il seguente significato: 

X: ascissa del punto noto sinistro deirintervallo di interpolazione, valore 
S 5 inizialmente corrispondente ad A; 

Y: ascissa del punto noto destro utilizzato per Tinterpoiazione, valore 
inizialmente corrispondente a B; 

Z: ascissa di un punto intermedio fra X e Y; 

W: ascissa del punto estremo sinistro, utilizzato effettivamente per 
20 Tinterpoiazione, corrispondente di volta in volta al punto X o al punto X-l ; 
f(k): ordinata corrispondente ad una generica ascissa k. 

Le variabili X e Y vengono inizializzate ponendole uguali agli estremi A e B, nel 
blocco iniziale indicato dal riferimento 32 in figura 3. 

Successivamente, nel blocco 36 in figura, si procede calcolando l'ascissa Z, 
25 intesa ovviamente come numero intero, del punto intermedio fra X e Y. H valore 
dell'ascissa Z viene arrotondato airintero inferiore, mediante la funzione FLOOR 
(indicata simbolicamente in figura). Nello stesso blocco si procede calcolando l'ascissa 
del punto estremo sinistro W, che risulta essere uguale a X o a X-l a seconda che il 
valore di X+Y sia rispettivamente pari o dispari (applicando la formula W=2«Z-Y). Si 
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calcola poi 1'ordinata f(Z) del punto intermedio Z facendo la media aritmetica fra 
1'ordinata del punto estremo sinistra W e 1'ordinata dell'estremo destro Y, che sono note. 

Si procede poi verificando, sempre all'interno del blocco 36, se il punto appena 
calcolato, di ascissa Z, corrisponde al punto di ascissa X+l. Se tale condizione non & 
5 verificata si assegna il valore di Z all'ascissa Y del punto noto destro, e si applicano 
nuovamente le operazioni contenute nel blocco 36. In pratica l'algoritmo del blocco 36 
viene applicato ricorsivamente sul semi-interval lo alia sinistra del punto 
precedentemente calcolato fino aquando non viene raggiunto il punto di ascissa X+l. 

Successivamente, vedi le operazioni indicate complessivamente con il 
10 riferimento 38, si procede per ascisse crescenti, partendo dairultimo punto calcolato 
X+l, cercando un primo punto ancora da calcolare. Una volta trovato tale punto, si 
procede eseguendo le operazioni contenute nel blocco 40, oppure si termina l'algoritmo 
se si e gih raggiunto l'estremo destro B dell'intervallo di interpolazione, poich<5 tutti i 
punti sono gii stati calcolati. 

15 Se invece b stato trovato uri punto ancora da calcolare, la cui ascissa e pari al 

valore assunto dalla variabile Z, si proceed .: ^v^.sado alia variabile X il valore Z-l, 
che rappresenta il primo punto noto a sinl^r cercando, per successivi tentativi, il v 
primo punto noto a destra di Z, valore che viene assegnato alia variabile Y. Per fare 
quest'ultima operazione la variabile Y viene incrementata di un'unitk per volta, fino a 

20 quando non si raggiunge un corrispondente coefficiente noto. Tali operazioni sono 
svolte airinterno del blocco 40. 

Con i nuovi valori di X e Y si applica nuovamente Talgoritmo ricorsivo del 
blocco 36, ed in cascata la procedura, contenuta nel blocco 38, di ricerca di un nuovo 
punto ancora da calcolare. Se tutti i punti sono stati calcolati, I'algoritmo termina, 
25 altrimenti si determina un nuovo intervallo X-Y a cui applicare nuovamente I'algoritmo 
ricorsivo 36. 

L'algoritmo precedentemente descritto, illustrato in figura 3, consente di 
calcolare, effettuando esclusivamente somme, divisioni per due ed operazioni di 
confronto fra registri, un qualsiasi numero di coefficient! di canale contenuti in un 
30 intervallo delimitato a sinistra da due coefficient noti (di ascissa A e A-l) e a destra da 
un coefficiente noto (di ascissa B). 
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Per meglio comprendere l'algoritmo di interpolazione generale precedentemente 
descritto, esaminiamo ora la sua applicazione a sei casi particolari, tipici del sistema 
UMTS, dal caso in cui i simboli pilota (Npilot) sono tre, ed i coefficient! da calcolare 
sono sette, fino al caso in cui i simboli pilota sono otto ed i coefficienti di canale da 
5 calcolare sono solamente due. 

II primo caso (N P ilot=3) e illustrate schematicamente in figura 4. Sull'ascissa 
sono indicati i coefficienti di canale corrispondenti ai tre simboli pilota dello slot 
corrente, rispettivamente alio posizioni di ascissa 0, 1 (A-l) e 2 (A), ed al primo 
simbolo pilota dello slot successivo, alia posizione di ascissa 10 (B). I coefficienti di 
1 0 canale da calcolare sono pertanto quelli corrispondenti alle posizioni di ascissa da 3 a 9. 

L'operazione di interpolazione avviene per passi successivi; il numero di passi 
dipende direttamente dal numero di coefficienti di canale da calcolare. Nel caso di 
figura 4 il calcolo viene completato in sette passi (Step 1 - Step 7). 

Al primo passo, Step 1, si applica una prima volta l'algoritmo del blocco 36 del 
1 5 diagramma di flussc m %ura 3 assegnando alle variabili X e Y rispettivamente i valori 
•di A (che vale 2.) <>.■• & &%vale 10): ■ • • ' . 

X=A=2 
Y=B=10 
Pertanto risulta: 

20 Z=FLOOR [(2+10)/2]=6 ; ascissa del punto intermedio da calcolare (indicato in 

figura con la lettera C); 

W=2-6-l0=2 ; ascissa del punto estremo sinistra, in questo caso corrispondente a 

X. 

Una volta note le ascisse del punto estremo sinistra W, del punto estremo destro 
25 Y e del punto intermedio Z di cui si vuole calcolare l'ordinata, si precede a calcolare il 
valore del corrispondente coefficiente C(6)=d(6)+jCQ(6), calcolando le component! C, 
e C Q con la regola della media aritmetica (punto C nel grafico di figura 4): 
C,(6)=[C,(2)+C,(lO)]/2 e C q (6)=[Cq(2)+Cq(10)]/2 
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Dopo aver verificato che l'ascissa Z del punto calcolato, in questo caso pari a 6, 
non corrisponde ancora al punto X+l=3, si assegna il valore Z=6 alia variabile Y e si 
calcola un nuovo coefficiente (punto D, step 2): 

X=2; Y=6 

5 Z=FLOOR [(2+6)/2]=4 

W=24-6=2 

Ci(4)=[Ci(2)+Ci(6)]/2 e C q (4)=[Cq(2)+C q (6)]/2 

Anche in questo caso la condizione di uscita dal blocco 36 non e verificata, in 
quanto Z=4 e di verso da X+l=3, pertanto si procede calcolando un nuovo coefficiente 
10 (punto E, step 3): 

X=2; Y=Z=4 

Z=FLOOR [(244)/2]=3 

W=2-3-4=2 

C,(3)=[C I (2)+C,(4)]/2 e C q (3)=[Cq(2)+Cq(4)]/2 

• > A questo punto si verifica la condizione Z=X+1 che consente di passare dal 

blocco 36 al blocco 38. 

Applicando le regole descritte nel blocco 38 trovo che il primo punto ancora da 
calcolare alia destra del punto di ascissa Z=3 e il punto di ascissa 5; verifico che tale 
punto non corrisponde all'estremo destro B=10 e proseguo nel blocco 40 avendo 
20 assegnato alia variabile Z il valore 5. 

Nel blocco 40 determino il primo punto noto a sinistra X=Z-l=4 ed il primo 
punto noto a destra Y=6, e con questi valori di X e Y ritorno al blocco 36 dove calcolo 
un nuovo coefficiente (punto F, step 4): 

X=4; Y=6 

25 Z=FLOOR [(4+6)/2]=5 

W=2-5-6=4 

Ci(5)=[Cr(4)+Q(6)]/2 e Cq(5)=[C q (4)+C q (6)]/2 

In questo caso la condizione Z=X+1 b subito verificata e si passa nuovamente al 
blocco 38 in cui si determina che il primo punto sulla destra ancora da calcolare & il 
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punto di ascissa Z=7. Nel blocco 40 determino il primo punto nolo a sinistra X=Z-1=6 
ed il primo punto noto a destra Y=10, e con questi valori di X e Y ritorno al blocco 36 
dove calcolo un nuovo coefficiente (punto G, step 5): 
X=6 ; Y=10 

Z=FLOOR [(6+10)/2]=8 
W=2-8-10=6 

C,(8)=[C,(6)+Ci(10)]/2 e C q (8)=[Cq(6)+Cq(10)]/2 

La condizione Z=X+1 non b verificata, pertanto resto all'interno del blocco 36 
per calcolare un nuovo coefficiente (punto H, step 6): 
X=6; Y=Z=8 
Z=FLOOR [(6+8)/2]=7 

W=2-7-8=6 

C,(7)=[C 1 (6)+Q(8)]/2 e Cq(7)=[C q (6)+C q (8)]/2 

Avendo raggiunto il punto Z=X+1=7 passo al blocco 38 in cui verv-^iS? 
aver ancora raggiunto l'estremo destroy ed in cui determino l'ascissa Z=9 ; M' A* v^d 
punto da calcolare. Nel blocco 40 determino il primo punto noto a sinistra X=^l=8 ed 
il primo punto noto a destra Y=10, e con questi valori di X e Y ritorno al blocco 36 
dove calcolo un nuovo coefficiente (punto I, step 7): 

X=8;Y=10 

Z=FLOOR [(8+10)/2]=9 
W=29-10=8 

QPMCW+COO)]^ e Cq(9)=[Cq(8)+Cq(10)]/2 

La condizione Z=X+1 risulta ora verificata, pertanto passo al blocco 38 in cui la 
variabile Z viene incrementata passando al valore Z=10 che, essendo corrispondente 
all'estremo destro B, mi porta a terminare Talgoritmo di interpolazione (passo STOP in 
ftgura 3). 

Come si pu6 notare dalle suddette equazioni tutti i coefficienti di canale sono 
stati calcolati effettuando solamente somme e divisioni per due, operazioni facilmente 
effettuabili direttamente in hardware. 
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Nel caso sopra descritto (N P1U yr=3) il coefficiente noto di ascissa A-l non 6 stato 
utilizzato, in quanto la distanza fra i punti A e B risulta essere una potenza di 2, in 
particolare 2 3 =8. Questa condizione si verifica anche nel caso N P iLo-r=7, in cui tale 
distanza 2 2 =4. 

L'algoritmo illustrato in figura 3 consente comunque di calcolare tutti i 
coefficienti intermedi anche nei casi in cui tale distanza non sia una poienza di due, 
come nei restanti casi qui di seguito illustrati. 

In figura 5 b illustrato il caso in cui i simboli pilota sono quattro (N Ptt oT=4) ed i 
coefficient! da calcolare sono sei, da C(4) fino a C(9). Lo schema di interpolazione 
illustrato in figura 5, come risultera dalla successiva spiegazione dettagliata, utilizza 
come estremo sinistro per 1'interpolazione il punto di ascissa A-l. Questa 
approssimazione, anche se introduce una Ieggera degradazione delle caratteristiche 
complessive del sistema, semplifica di molto il calcolo dei coefficienti interpolati. 

Vediamo ora come si applica l'algoritmo di figura 3 al presente caso (N P ilot=4). 

Al primo passo, step I, si applica una prim?.. v-D'ta. i'algoritmo del blocco 36 del * 
diagramma di flusso di figura 3 assegnando all© Ur^o Z g:.e Y rispetti vamehte i valori 
di A (che vale 3) e di B (che vale 10): 

X=A (A=3) 

Y=B (B=10) 

Pertanto risulta: 

Z=FLOOR [(3+10)/2]=6 ; ascissa del punto intermedio che si vuole calcolare 
(indicato in figura con la lettera C); 

W=2-6-10=2 ; ascissa del punto estremo sinistro, in questo caso non 
corrispondente a X. 

Una volta note le ascisse del punto estremo sinistro W, del punto estremo destro 
Y e del punto intermedio Z di cui si vuole calcolare l'ordinata, si precede a calcolare il 
valore del corrispondente coefficiente C(6)=Ci(6)+jCQ(6), calcolando le componenti Q 
c Cq con la regola della media aritmetica (punto C nel grafico di figura 5): 

Ci(6)=rC,(2)+C,(lO)]/2 e C Q (6)=[C Q (2)+C Q (10)]/2 
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Dopo aver veriftcato che l'ascissa Z del punto calcolato, in questo caso pari a 6, 
non corrisponde ancora al punto X + l=4, si assegna il valore Z=6 alia variabile Y e si 
calcola un nuovo coefficiente (punto D, step 2): 
X=3; Y=6 

Z=FLOOR [(3+6)/2]=4 
W=24-6=2 

C,(4)=[C,(2)+Ci(6)]/2 e C Q (4)=[C Q (2)+C Q (6)]/2 

In questo caso la condizione di uscita dal blocco 36 risulta verificata, in quanto 
Z=4 e pari a X+l=4, pertanto si passa al blocco 38. 

Applicando le regole descritte nel blocco 38 trovo che il primo punto ancora da 
calcolare alia destra del punto di ascissa 2=4 e il punto di ascissa 5; verifico che tale 
punto non corrisponde all'estremo destro B=10 e proseguo nel blocco 40 avendo 
assegnato alia variabile Z il valore 5. 

Nel blocco 40 determine il primo punto noto a sinistra X=Z-1=4 ed il primo 
punto noto a destra YvJjO-i <;,,,<, questi valori di X e Y ritorno al blocco 36 dove calcplo . v 
un nuovo coefficienie" j :«cB, step 3): .- ~- 

X=4; Y=6 

Z=FLOOR [(4+6)/2]=5 
W=2- 5-6=4 

20 C I (5)=[C I (4)+Q(6)]/2 e C Q (5)=[C Q (4)+C Q (6)]/2 

A questo punto si verifica la condizione Z=X + 1 che consente di passare dal 
blocco 36 al blocco 38. 

Applicando le regole deseritte nel blocco 38 trovo che il primo pnnto ancora da 
calcolare alia destra del pnnto di ascissa *=5 e il punto di ascissa 7; verifico che tale 
25 punto non corrisponde all'estremo detfro B=10 e prosegno nel blocco 40 avendo 
assegnato alia variabile Z il valore 7. 

Nel blocco 40 determino il primo punto noto a sinistra X=Z-1=6 ed il primo 
punto noto a destra Y=10, e con questi valori di X e Y ritorno al blocco 36 dove calcolo 
un nuovo coefficiente (punto F, step 4): 
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X=6; Y=IO 

Z=FLOOR [(6+10)/2]=8 
W=2-8-10=6 

Ci(8)=[Q(6)+C,(10)]/2 e C Q (8)=[C Q (6)+C Q (10)]/2 

5 La condizione Z=X+1 non e verificata, pertanto resto aH'interno del blocco 36 

per calcolare un nuovo coefficiente (punto G, step 5): 

X=6; Y=Z=8 

Z=FLOOR [(6+8)/2]=7 

W=2-7-8=6 

10 C,(7)=[C,(6)+C,(8)]/2 e Cq(7)=[C q (6)+C q (8)]/2 

In questo caso la condizione Z=X+1 e subito verificata e si passa nuovamente al 
blocco 38 in cui si determina che il prime- punto sulla destra ancora da calcolare e il 
punto di ascissa Z=9. Nel blocco 40 determino il primo punto noto a sinistra X=Z-1=8 
il primo punto noto a destra Y=10, e con questi valori di X e Y ritorno al blocco 36 
1 > : calcolo un riuovo coefficiente (punto H, step 6): 

X=8; Y=10 

Z=FLOOR [(8+10)/2]=9 
W=2-9-10=8 

Ci(9)=[Ci(8)+C,(10)]/2 e Cq(9)=[Cq(8)+Cq(10)]/2 

20 La condizione Z=X+l risulta ora verificata, pertanto passo al blocco 38 in cui la 

variabile Z viene incrementata passando al valore Z=10 che, essendo corrispondente 
all'estremo destro Y, mi porta a terminare ralgoritmo di interpolazione (passo STOP in 
figura 3). 

La figura 6 illustra il caso in cui ciascuno slot contiene cinque simboli pilota 
25 (N P ilot=5) ed il numero di coefficients di canale da calcolare e pari a cinque. 
L'operazione viene svolta utilizzando sia il coefficiente C(4) corrispondente airultimo 
simbolo pilota sia il coefficiente C(3) corrispondente al penultinio simbolo pilota. 
Anche in questo caso Tapprossimazione introduce una leggera degradazione delle 
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caratteristiche complessive del sistema, ma il calcolo dei coefficienti interpolati risulta 
essere di molto semplificato. 

I cinque passi necessari per il calcolo dei cinque coefficienti verranno ora 
descritti in forma abbreviate riportando i valori assunti dalle variabili, dato che la 
5 metodologia di applicazione dell'algoritmo di figura 3 risulta essere equivalente a quella 
descrittacon riferimento ai casi precedenti. 

Nel primo passo si calcola il coefficiente C(7) (punto C, step 1)3 

X =4; Y=10 

Z=FLOOR [(4+10)/2]=7 

10 W=2-7-10=4 

C,(7)=[C,(4)+Ci(10)]/2 e Cq(7)=[C q (4)+Cq(10)]/2 

Nel secondo passo si calcola il coefficiente C(5) (punto D, step 2): 

X=4; Y=7 

Z=FLOOR [(4+7)/2]=5 

15 W=2- 5-7=3 

Ci(5)=[C,(3)+C I (7)]/2 e C q (5)=[Cq(3)+C q (7)]/2 

Nel terzo passo si calcola il coefficiente C(6) (punto E, step 3): 

X=5; Y=7 

Z=FLOOR [(5+7)/2]=6 

20 W=2-6-7=5 

C,(6)=[C,(5)+C I (7)]/2 e C q (6)=[Cq(5)+C q (7)]/2 

Nel quarto passo si calcola il coefficiente C(8) (punto F, step 4): 

X=7;Y=10 

Z^FLOOR [(7+10)/2]=8 

25 W=2-8-10=6 

C,(8)=LC,(6)+C 1 (10)]/2 e C Q (8)=[C Q (6)+C Q (10)]/2 

Nel quinto passo si calcola il coefficiente C(9) (punto G, step 5): 

X=8 ; Y=10 
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Z=FLOOR [(8+10)/2]=9 
W=2-9-10=8 

Ci(9)=[C,C8)+C,(10)]/2 e Cq(9)=[Cq(8)+Cq(10)]/2 

Anche in questo caso i cinque coefficienti C(5) ... C(9) sono stati ottenuti 
effettuando esclusivamente somme e divisioni per due. 

H caso in cui ciascuno slot contiene sei simboli pilota (Npilot=6) ed il numero di 
coefficienti di canale da calcolare e pari a quattro e invece illustrato in figura 7. Anche 
in questo caso si utilizza sia il coefficiente C(5), corrispondente all*uItimo simbolo 
pilota, sia il coefficiente C(4), corrispondente al penultimo simbolo pilota. 

Lo schema di interpolazione di figura 7 e basato su quattro passi consecutivi, qui 
di seguito riassunti: 

Nel primo passo si calcola il coefficiente C(7) (punto C, step 1): 

X=5; Y=10 

Z=FLOOR [(5+10)/2]=7 
W=2-7-10=4 

C,(7)=[C,(4)+C,(10)]/2 e Cq(7)=[Cq(4)+Cq(10^7'2 

Nel secondo passo si calcola il coefficiente C(6) (punto D, step 2): 

X=5; Y=7 

Z=FLOOR [(5+7)/2]=6 
W=2-6-7=5 

Ci(6)=[Ci(5)+C,(7)]/2 e C q (6)=[C q (5)+Cq(7)]/2 

Nel terzo passo si calcola il coefficiente C(8) (punto E, step 3): 

X=7; Y=10 

Z=FLOOR [(7+10)/2]=8 
W=2-8-10=6 

Cl(8)=[Ci(6)+C,(10)]/2 e C Q (8)=[C Q (6)+C Q (10)]/2 

Nel quarto passo si calcola il coefficiente C(9) (punto F, step 4): 

X=8 ; Y=10 
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Z=FLOOR [(8+10)/2J=9 
W=2-9-l0=8 

C,(9)=[C,(8)+C,(10)]/2 e Cq(9)=[Cq(8)+Cq(10)]/2 

Anche in queslo caso le operazioni effettuate sono solamente somme e divisioni 
5 per due. 

In figura 8 e illustrato il caso in cui lo slot contiene sette simboli pilota 
(N PIU yr=7) ed il numero di coefficienti di canale da calcolare e pari a tre. L'operazione 
viene svolta in tre passi, utilizzando il coefficiente C(6) corrispondente all'ultimo 
simbolo pilota. 

10 Ncl primo passo si calcola il coefficiente C(8) (punto C, step 1): 

X=6; Y=10 

Z=FLOOR [(6+10)/2]=8 
W=2-8-10=6 

C^CC^+CCIO)],^ • >l Cq(8)=[Cq(6)+Cq(10)]/2 - 
15 Nel secondo pas? " .'. ^ coefficiente G(7) (punto D, step 2): 

X=6; Y=8 

Z=FLOOR [(6+8)/2]=7 
W=2- 7-8=6 

C,(7)=[C,(6)+C,(8)]/2 e Cq(7)=[Cq(6)+Cq(8)]/2 
20 Nel terzo passo si calcola il coefficiente C(9) (punto E, step 3): 

X=8 ; Y=10 

Z=FLOOR [(8+10)/2]=9 
W=2-9-10=8 

C,(9)=[C,(8)+C,(10)]/2 e C Q (9)=[C Q (8)+C Q (10)]/2 
25 In figura 9 b illustrato invece 1'ultimo caso, in cui ciascuno slot contiene otto 

simboli pilota (N p ,lot=8) ed il numero di coefficienti di canale da calcolare e pari a due. 
Questa operazione viene svolta in due passi, utilizzando il coefficiente C(6) che 
corrisponde, in questo caso, al penultimo simbolo pilota. 



16 



Nel primo passo si calcola il coefficiente C(8) (punto C, step 1): 
X=6; Y=10 

Z=FLOOR [(6+10)/2]=8 
W=2-8-10=6 

C|(8)=[Q(6)+Ci(10)]/2 e C q (8)=[Cq(6)+C q (10)]/2 

Nel secondo passo si calcola il coefficiente C(9) (punto D, step 2): 

X=8; Y=10 

Z=FLOOR [(8+10)/2]=9 
W=2-9-10=8 

Ci(9MCi(8)+Q(10)]>2 e Cq(9)=[Cq(8)+C q (10)]/2 

Anche in questo caso il fatto di utilizzare il coefficiente di canale corrispondente 
al penultimo simbolo pilota introduce una leggera degradazione delle caratteristiche 
complessive del sistema, ma il calcolo dei coefficienti interpolati risulta essere di molto 

I;>ei casi precedentemente descritti ed illustrati con riferimento alle figure da 4 a 
9 possono essere schematizzati nel diagramma a stati illustrato in figura 10. Tl 
diagramma di figura 10, la cui sequenza di operazioni dipende essenzialmente dal 
valore del parametro Npilot, illustra come sia possibile implementare, con una semplice 
macchina a stati, un dispositivo hardware o software che realizza il metodo della 
presente invenzione. Nel diagramma b rappresentata solo la componente in fase Ci(k) 
dei coefficienti di canale, in quanto le formule per il calcolo della corrispondente 
componente in quadratura Cq(1c) risultano essere equivalent! 

In figura 11 e invece illustrato un esempio di applicazione, riferito al caso 
Npilot=6, di un secondo metodo di interpolazione realizzato secondo la presente 
invenzione. Nell'esempio di figura 11 per calcolare una plurality di coefficienti di 
canale, corrispondenti alle posizioni di ascissa da 6 a 9, vengono utilizzati i coefficienti 
di canale noti corrispondenti rispettivamente aH'ultimo simbolo pilota, di ascissa A, 
dello slot corrente L ed ai primi due simboli pilota, di ascisse B e B+l, dello slot 
successivo Lr*-1 a quel lo corrente. 
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Per implementare questo secondo metodo di interpolazione e necessario 
memorizzare un simbolo pilota in piu rispetto ai casi precedentemente descritti, quello 
di ascissa 1 1 (B+l), corrispondente al secondo simbolo pilota dello slot successivo. 

Per comprendere la sequenza di passi illustrata nell'esempio di applicazione di 
5 figura 1 1, e per capire come lo stesso metodo possa essere applicato anche ai restanti 
casi, da N Pnj oT=3 a N P ilot=8, analizziamo ora in dettaglio il secondo metodo di 
interpolazione realizzato secondo la presente invenzione. 

Questo secondo metodo consente di calcolare per interpolazione una pluralita di 
punti compresi fra un ultimo coefficiente di canale, di ascissa A, di uno slot corrente L, 
10 ed i primi due coefficienti di canale, di ascisse B e B+ 1 , di uno slot successivo L+ 1 . 

Supponiamo di utilizzare, all*interno di questo algoritmo alcune variabili che 
hanno lo stesso significato delle variabili precedentemente utilizzate con riferimento 
all'algoritmo di figura 3: 

X: ascissa del punto noto sinistra utilizzato per I'interpolazione, valore 
15 inizialmente corrispondente ad A; 

Y: ascissa del punto noto destro dell'intervallo di interpolazione, va3*Mst; 
inizialmente corrispondente a B; 

Z: ascissa di un punto intermedio fra X e Y; 

W: ascissa del punto estremo destro, utilizzato effettivamente per 
20 I'interpolazione, corrispondente di volta in volta al punto Y o al punto Y+ 1 ; 

f(k): ordinata corrispondente ad una generica ascissa k. 

Le variabili X e Y vengono inizializzate ponendole uguali agli estremi A e B. 

Successivamente si precede calcolando l'ascissa Z, intesa ovviamente come 
numero intero, del punto intermedio fra XeY.D valore dell'ascissa Z viene arrotondato 
25 all'intero superiore, mediante la funzione CEIL [(X+Y)/2]. Si procede poi calcolando 
l'ascissa del punto estremo destro W, che risulta essere uguale a Y o a Y+l a seconda 
che il valore di X+Y sia rispettivamente pari o dispari, utilizzando la formula W=2»Z- 
X. Si calcola poi l'ordinata f(Z) del punto intermedio Z facendo la media aritmetica fra 
l'ordinata del punto estremo destro W e l'ordinata dell'estremo sinistra X, che sono note. 
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Si precede poi verificando se il punto appena calcolato, di ascissa Z, corrisponde 
al punto di ascissa Y-l. Se tale condizione non e verificata si assegna il valore di Z alia 
variabile X del punto noto sinistra, e si applicano nuovamente le operazioni del 
paragrafo precedente. In pratica le suddette operazioni vengono applicate 

5 ricorsivamente su ciascun semi-intervallo alia destra del punto precedentemente 
calcolato fmo a quando non viene raggiunto il punto di ascissa Y-l . 

Successivamente si procede per ascisse decrescenti, partendo dalMtimo punto 
calcolato Y-l e cercando un primo punto ancora da calcolare. Se si raggiunge I'estremo 
sinistro A dell'intervallo si termina l'algoritmo, poiche tutti i punti sono gia stati 

10 calcolati. 

Se invece b stato trovato un punto ancora da calcolare, la cui ascissa e pari al 
valore assunto dalla variabile Z, si procede assegnando alia variabile Y il valore Z+l, 
che rappresenta il primo punto noto a destra, e si cerca, per successivi tentativi, il primo 
punto noto a sinistra di Z, valore che viene assegnato alia variabile X. Per fare 
1 5 quest'ultima operazione la variabile X viene decrementata di un'unita per volta, fino a - 
quando non si raggiunge un corrispondente coefficiente n*:t> 

Con i nuovi valori di X e Y si applica nuovamente i^Syittrnd ricorsivo per 41 
calcolo di un punto intermedio di ascissa Z, ed in cascata la procedura di ricerca di un 
nuovo punto ancora da calcolare. Se tutti i punti sono stati calcolati, l'algoritmo termina, 
20 altrimenti si determina un nuovo intervallo X-Y a cui applicare il suddetto algoritmo 
ricorsivo. 

Analizziamo ora come viene applicato l'algoritmo precedentemente descritto alio 
schema di inlerpolazione di figura 11. 

Nel primo passo si calcola il coefficiente C(8) (punto C, step 1): 
25 X=A=5; Y=B=10 

Z=CEIL [(5+10)/2]=8 
W=2-8-5=l 1 

Ci(8)=[C,(5)+C,(ll)]/2 e Cq(8)=[C q (5)+Cq(1 l)]/2 
Nel secondo passo si calcola il coefficiente C(9) (punto D, step 2): 
30 X=8;Y=10 
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Z=CE1L [(8+10)/2]=9 
W=29-8=10 

C,(9)=[C,(8)+C,(10)]/2 e Cq(9)=[Cq(8)+Cq(10)]/2 

Nel Lerzo passo si calcola il coefficiente C(7) (punto E, step 3): 

X=5; Y=8 

Z=CEIL [(5+8)/2]=7 . 

W=2-7-5=9 

C I (7)=[C,(5)+C,(9)]/2 e C Q (7)=[C Q (5)+C Q (9)]/2 

Nel quarto passo si calcola il coefficiente C(6) (punto F, step 4): 

X=5; Y=7 

Z=FLOOR [(5+7)/2]=6 
W=2-6-5=7 

C,(6)=[Ci(5)+Ci(7)]/2 e C^/^)=[C^(5)+Cq(7)]/2 

L'algoritmo descritto l 4 ■ ?>3f «* tf^tnetodo di interpolazione pu6 essere yisto 
corne una versione "speculare" deiraigoritmo precedentemente descritto con riferimento 
alia figura 3. Anch'esso infatti consente di calcolare, effettuando esclusivamente 
somme, divisioni per due ed operazioni di confronto fra registri, un qualsiasi numero di 
coefficienli di canale contenuti in un intervallo delimitato a sinistra da un coefficiente 
noto (di ascissa A) ed a destra da due coefficienti noti (di ascisse B e B+l). 

I due metodi di interpolazione precedentemente descritti possono essere 
utilizzati, in combinazione, per implementare un terzo metodo, che pud risultare 
vantaggioso in termini di velocity di esecuzione. 

Nel caso in cui siano disponibili, come estremi noti deirintervallo di 
interpolazione, due coefficienti di canale alia sinistra, di ascissa A-l e A, e due 
coefficienti di canale noti alia destra, di ascissa B e B+l, risulta possibile applicare "in 
parallelo" entrambi i metodi di interpolazione precedentemente descritti, rispettivamente 
in un sottointervallo sinistra ed in un sottointervallo destro dell'intervallo originale. E' 
sufficiente infatti, come primo passo, calcolare un punto medio C fra i due estremi A e 
B deirintervallo, arrotondando eventualmente l'ascissa di tale punto C all'intero 
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inferiore o all'intero superior©, indifferentemente, per definire un sottointervallo sinistro, 
compreso fra i punti di ascissa A-l, A e C, ed un sottointervallo destro compreso fra i 
punti di ascissa C e B, B+l. A questo punto e possibile applicare in parallelo, sfruttando 
quindi meglio le risorse hardware impegnate, il primo metodo di interpolazione nel 
5 sottointervallo sinistro ed il secondo metodo di interpolazione nel sottointervallo destro, 
giungendo quindi parallelamente al calcolo di tutti i coefficienti compresi nell'intervallo. 

Analogamente b possibile prevedere un'ulteriore variante dei metodi realizzati 
secondo la presente invenzione; ad esempio, e possibile utilizzare per l'interpolazione, 
anzich6 il coefficiente di canale noto immediatamente precedente all'ultimo coefficiente 
10 nolo dello slot corrente o immediatamente successivo al primo coefficiente noto dello 
slot successivo, altri coefficienti noti piu distanti dall'intervallo di interpolazione o una 
combinazione lineare degli stessi. 

Come esempio d'utilizzo di una combinazione lineare di coefficienti noti si 
considera il caso di tre piloti (N P iuyr=3). Dapprima si esegue una combinazione lineari 
15 dei tre piloti noti del gruppo di sinistra 
•.„,-.C I {2)/2+[Ci(0)+ C,(I)]/4 
CLq=Cq(2)/2+[C q (0)+ Cq(1)]/4 
e si assegna il risultato al coefficiente C(2). 
Q(2)=CLi 
20 Cq(2)= CLq 

II coefficiente C(l) viene anch'esso calcolato come combinazione lineare dei 
medesimi piloti: 

CL,=C,(l)/2+[C,(0)+ Ci(2)]/4 
CLq=Cq(1)/2+[C q (0)+ Cq(2)]/4 
25 II risultato si assegna al coefficiente C( 1). 

C,(1)=CL, 
Cq(1)=CLq 

Una combinazione lineare simile pub essere calcolata sulla base dei piloti dello 
slot successivo. 
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CLi=Ci(10)/2+[Q(1 1)+ C,(12)]/4 
CLq=C q (10)/2+[C q (1 1)+ C q (12)]/4 
e quindi assegnata al coefficiente C(10) 
Q(10)= CL, 
5 C 0 (10)=CL Q 

Le combinazioni lineari sono scelte in modo da essere implementate come 
somme e shift di bit analogamente a quanto illustrate per l'interpolazione. 

L'algoritmo di interpolazione per il calcolo dei coefficienti mancanti, pu6 quindi 
essere applicato seguendo lo schema gia illustrato per tre piloti e riportato in Figura 4. 
10 Questa variante del metodo trova la sua giustificazione nel caso molto frequente che i 
piloti siano affetti da rumore: una combinazione lineare, che rappresenta una media 
pesata dei coefficienti di canale noti ma affetti da rumore, fornisce una stima piil 
affidabile dei punti noti per l'interpolazione. 

La figura 12 illustra, in uno schema a blocchi, una possibile implementazione 
1 5 hardware di un dispositivo di stima di canale realizzato secondo la presente irivenziomv 
nel caso di una tratta di trasmissione verso una stazione radio base ed uno "spreading 
factor" del canale DPDCH pari a 256. 

II dispositivo di figura 12, definito Unita di Stima del Canale (CEU, Channel 
Estimation Unit), consente di estendere la stima del canale ai coefficienti di canale 
20 successivi ai simboli pilota, attuando uno dei metodi di interpolazione precedentemente 
descritti. 

L'unitk CEU riceve in ingresso i simboli di canale ricevuti dopo aver effettuato 
le operazioni di "descrambling", "despreading" con il codice di canalizzazione DPCCH 
e integrazione. I simboli pilota ricevuti che alimentano l'unita CEU sono indicati con Zi 
25 (k), dove il pedice "i u indica la componente di segnale (i = I per la componente in fase o 
i = Q per la componente in quadratura). 

La prima operazione fatta dall'unita. CEU e la moltiplicazione dei simboli pilota 
ricevuti per i simboli pilota di riferimento Xq (k), i simboli pilota di riferimento sono 
infatti noti al ricevitore. Successivamente b necessario invertire il segno della 
30 componente in fase. L'inversione di segno della componente in fase viene effettuata 
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moltiplicando per il fattore Mi, ponendo Mi = -1 per la componente in fase e Mi = +1 
per la componente in quadratura, come mostrato nello schema di figura 12. 

Una volta effettuate le suddette operazioni l'effetto dell a sequenza dei simboli 
pilota b slato rimosso, ed i valori complessi ottenuti rappresentano solamente lo 
5 spostamento di fase introdotto dal canale. Questi valori vengono memorizzati nella 
memoria 100 attraverso Tingresso contrassegnato in figura 10 come Input_port_2. 

Quando tutti i simboli pilota dello slot corrente ed uno o due simboli pilota 
iniziali dello slot successivo, a seconda del metodo che si vuole implernentare, sono 
stati ricevuti e memorizzati nella memoria 100, inizia il processo di interpolazione sotto 
10 il controllo dell'unite logica di controllo 102, (CLU, Control Logic Unit). 

L'unifo logica di controllo 102 b una macchina a stati finiti (FSM, Finite State 
Machine), che esegue una sequenza predefinita di operazioni, dipendente dal parametro 
di ingresso Npilot (corrispondente al numero di simboli pilota presenti nello slot 
corrente). Nel caso ad esempio del metodo di interpolazione precedentemente illustrate 
15 con riferimento alia figura 3, I'unilSt logica di controllo 102 z\?r\t*. la sequenza Ml? 
operazioni illustrata nel diagramma di stato di figura 10. 

L'operazione di base effettuata dalFunitfr logica di coiit-^o 102 e divisa nei r 
seguenti tre passi: 

- leggere il primo operando dalla memoria 100 e caricarlo nel primo registro 104; 

20 - leggere il secondo operando dalla memoria 100 e caricarlo nel secondo registro 106; 

- scrivere la media aritmetica dei due operandi nella memoria 100 attraverso il suo 
ingresso Input_jport_l . 

II calcolo della media aritmetica dei due operandi richiede una somma e una 
divisione per due, la somma viene effettuata nel sommatore 108 e la divisione per due 
25 nel blocco 110 che effettua uno shift a destra, in pratica scartando il bit meno 
significativo (LSB, Less Significant Bit) del dato uscente dal sommatore 108. 

Tutti i calcoli eseguiti richiedono solo somme e shift a destra, operazioni che 
sono facilmente implementabili in hardware. 

Quando tutti i coefficienti di canale Q(k), con k che va da 0 a 9, sono stati 
30 calcolati, questi vengono sequenzialmente letti dalla memoria e forniti all'unife di 
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compensazione di canaJe (CCU, Channel Comensation Unit), non illustrata in figura 10 
in quanto la sua funzione b ben nota al tecnico del ramo. 

II fattore moltiplicativo indicato in figura 12 come M 2 6 necessario per invertire 
il segno della componente in quadratura e per fornire il complesso coniugato delle stirae 
di canale all'unita di compensazione di canale (CCU). M2 e pari a +1 per la componente 
in fase ed e pari a -1 per la componente in quadratura. 
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RTVENDICAZIONI 

1. Metodo per la stima della fiinzione di trasferimento di almeno un canale di 
Irasmissione in un sistema di ricezione per reti di comunicazione, mediante il calcolo, 
per mezzo di un algoritmo di interpolazione, di una plurality di coefficienti di canale 
compresi fra due gruppi di coefficienti di canale noti, ciascun coefficiente di canale 
essendo associato ad un valore di ascissa intero su un asse temporale, caratterizzato dal 
fatto che il calcolo di detta plurality di coefficienti di canale viene effettuato applicando 
ripetutamente un algoritmo di interpolazione atto a calcolare un punto intermedio (Z, 
f(Z)) fra un primo estremo ed un secondo estremo di un determinato intervallo, detto 
primo estremo essendo costituito da almeno due punti noti e detto secondo estremo 
essendo costituito da almeno un punto nolo, detto punto intermedio avendo per ascissa 
(Z) il valore di ascissa del punto medio fra i punti che delimitano detto intervallo 
arrotondato airintero pit vicino a detto primo estremo, ed avendo per ordinata (f(Z)) la 
media aritmetica fra l'ordinata del punto noto di detto secondo estremo e l'ordinata di un 
punto, scelto fra i due punti noti di detto primo estremo, avente una di stanza da detto 
punto *r£> o- ^ su detto asse temporale, pari alia distanza fra detto punto intermedio 
ed il pu.i r. . ^i<ietto secondo estremo. > 

2. Metodo secondo la rivendicazione 1, in cui detti coefficienti di canale da 
calcolare sono compresi fra un primo coefficiente di canale noto, di ascissa A, 
corrispondente ad un ultimo simbolo pilota di uno slot corrente (L) ed un secondo 
coefficiente di canale noto, di ascissa B, corrispondente ad un primo simbolo pilota di 
uno slot (L+l) successivo a detto slot corrente, essendo inoltre noto un terzo 
coefficiente di canale, di ascissa A-l, alia sinistra di detto primo coefficiente di canale 
di ascissa A, ed il calcolo di detti coefficienti di canale avviene mediante i seguenti 
passi: 

a) applicare ripetutamente in modo ricorsivo detto algoritmo di interpolazione 
nell'intervallo delimitato da detti coefficienti di canale noti di ascissa A e B, effettuando 
una prima iterazione in cui viene calcolato un primo coefficiente intermedio (di ascissa 
C), ed effettuando successive iterazioni dello stesso algoritmo in sotto-intervalli di volta 
in volta delimitati a sinistra da detto coefficiente di canale noto di ascissa A ed a destra 
dal coefficiente intermedio calcolato nella precedente iterazione, fino a quando non e 
stato raggiunto e calcolato il punto di ascissa A+l ; 
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b) cercare, procedendo per ascisse crescenti, un primo punto ancora da calcolare alia 
destra dell'ultimo coefficiente intermedio calcolato, definire come estremi di un nuovo 
intervallo di applicazione di detto algoritmo di interpolazione rispetlivamente il prirno 
punto noto a sinistra ed il primo punto noto a destra rispetto a detto punto ancora da 

5 calcolare, ed applicare nuovamente, in modo ricorsivo, detto algoritmo di interpolazione 
su detto nuovo intervallo effettuando successive iterazioni dello stesso algoritmo in 
sotto-intervalli di volta in volta delimitati dairestremo sinistra di detto nuovo intervallo 
e da un cstremo destro costituito dal coefficiente intermedio calcolato nella precedente 
iterazione, fino a quando non e stato raggiunto e calcolato il punto immediatamente 
10 adiacente all'estremo sinistra di detto nuovo intervallo; 

c) ripetere il passo b) fino a quando non b stato calcolato il coefficiente di canale 
associato al valore di ascissaB-1. 

3. Metodo secondo la rivendicazione 2, in cui ciascuno slot contiene tre simboli 
pilota (0, 1, 2), detto primo coefficiente di canale noto di ascissa A & il coefficiente 

15 C(2)=Ci(2)+Cq(2) corrispondente all'ultimo simbolo pilota (2) dello slot corrente (L), 
detto secondo coefficiente di canale noto di ascissa B & il coefficiente 
C(10)=Ci(10)+Cq(10) corrispondente al primo simbolo pilota (10) di uno slot 
successivo (L+l), e detto terzo coefficiente di canale noto di ascissa A-l £ il 
coefficiente C(1)=Ci(1)+Cq(1) corrispondente al penultimo simbolo pilota (1) dello slot 

20 corrente (L), ed il calcolo dei coefficienti di canale C(k)=d(k)+CQ(k), con k=3..9, 
avviene secondo la sequenza: 

C,(6)=[Ci(2)+C,(10)]/2 ; C q (6)=[C q (2)+Cq(10)]/2 ; 

C,(4)-[C I (2)+C,(6)]/2 ; Cq(4)=[Cq(2)+C q (6)1/2 ; 

Q(3)=[Ci(2)+Q(4)]/2 ; C Q (3)=[C Q (2)+C Q (4)]/2 ; 

25 Ci(5)=[Q(4)+C,(6)]/2 ; C Q (5)=[C Q (4)+C Q (6)]/2 ; 

Q(8)=[C I (6)+Ci(10)]/2 ; Cq(8)=[C q (6)+Cq(10)]/2 ; 

C,(7)=[C,(6)+C,(8)]/2 ; Cq(7)=[Cq(6)+Cq(8)]/2 ; 

C,(9)=[Ci(8)+Q(10)]/2 ; Cq(9)=[Cq(8)+C q (I0)]/2 . 

4. Metodo secondo la rivendicazione 2, in cui ciascuno slot contiene quattro 
30 simboli pilota (0, 1, 2, 3), detto primo coefficiente di canale noto di ascissa A e il 

coefficiente C(3)=Q(3)+Cq(3) corrispondente airultimo simbolo pilota (3) dello slot 
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correnle (L), detto secondo coefficiente di canale noto di ascissa B e il coefficiente 
C(10)=Q(10)+Cq(10) corrispondente al primo simbolo pilota (10) di uno slot 
successi vo (L+l), e detto terzo coefficiente di canale noto di ascissa A-l e il 
coefficiente C(2)=Ci(2)+Cq(2) corrispondente al penultimo simbolo pilota (2) dello slot 
corrente (L), ed il calcolo dei coefficienti di canale C(k)=Ci(k)+C Q (k), con k=4..9, 
avviene secondo la sequenza: 

C,(6)=[C,(2)+C I (10)]/2 ; Cq(6)=[Cq(2)+Cq(10)]/2 ; 

C,(4)=[C,(2)+C 1 (6)]/2 ; Cq(4)=[C q (2)+Cq(6)]/2 ; 
C,(5)=[C,(4)+C,(6)]/2 ; Cq(5)=[C q (4)+Cq(6)]/2 ; 
C,(8)=[C,(6)+Ci(10)]/2 ; Cq(8)=[C q (6)+C q (10)]/2 ; 
Ci(7)=[Ci(6)+C,(8)]/2 ; Cq(7)=[C q (6)+C q (8)]/2 ; 
C,(9)=[d(8)+C,(10)]/2 ; Cq(9)=[C q (8)+Cq(10)]/2 . 

5. Metodo secondo la rivendicazione 2, in cui ciascuno slot contiene cinque 
simboli pilota (0, 1, 2, 3, 4), detto primo coefficiente di canale noto di ascissa A e il 
coefficiente C(4)=Q(4)+Cq(4) corrispondente all'uUimo simbolo V} (4? dello slot 
corrente (L), detto secondo coefficiente di canale noto di ascL ; - '■:> ? : w^Oefficiente 
C(10)=C 1 (10)+Cq(10) corrispondente al primo simbolo pilota (10) di uno slot 
successivo (L+l), e detto terzo coefficiente di canale noto di ascissa A-l e il 
coefficiente C(3)=Q(3)+C Q (3) corrispondente al penultimo simbolo pilota (3) dello slot 
corrente (L), ed il calcolo dei coefficienti di canale C(k)=Ci(k)+C Q (k), con k=5..9, 
avviene secondo la sequenza: 

C,(7)=[C,(4)+C,(10)]/2 ; C q (7)=[C q (4)+Cq(10)]/2 ; 

Ci(5)=rC,(3)+C,(7)]/2 ; Cq(5)=[Cq(3)+Cq(7)]/2 ; 

C,(6)=[C,(5)+Ci(7)]/2 ; Cq(6)=[Cq(5)+C q (7)]/2 ; 

C,(8)=[C,(6)+C,(10)]/2 ; C q (8)=[Cq(6)+Cq(10)]/2 ; 

C i (9)=[C,(8)+C,(10)]/2 ; Cq(9)=[C q (8)+C q (10)]/2 . 

6. Metodo secondo la rivendicazione 2, in cui ciascuno slot contiene sei simboli 
pilota (0, 1, 2, 3, 4, 5), detto primo coefficiente di canale noto di ascissa A e il 
coefficiente C(5)=Ci(5)+C Q (5) corrispondente all'ultimo simbolo pilota (5) dello slot 
corrente (L), detto secondo coefficiente di canale noto di ascissa B 6 il coefficiente 
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C(IO)=Ci(10)+Cq(10) corrispondente al primo simbolo pilota (10) di uno slot 
successi vo (L+l), e detto terzo coefficiente di canale noto di ascissa A-l e il 
coefficiente C(4)=Q(4)+Cq(4) corrispondente al penultimo simbolo pilota (4) dello slot 
corrente (L), ed il calcolo dei coefficient di canale C(k)=Ci(k)+C Q (k), con k=6„9, 
5 avviene secondo la sequenza: 

C,(7)=[C,(4)+C,(10)]/2 ; Cq(7)=[Cq(4)+Cq(10)]/2 ; 

C,(6)=[d(5)+C,(7)]/2 ; C q (6)=[C q (5)+Cq(7)]/2 ; 

C ] (8)=[C,(6)+C,(10)]/2 ; C Q (8)=[C Q (6)+C Q (10)]/2 ; 

C I (9)=[C l (8)+C,(10)]/2 ; C q (9)=[Cq(8)+Cq(10)]/2 . 

10 7. Metodo secondo la rivendicazione 2, in cui ciascuno slot contiene sette simboli 
pilota (0, 1, 2, 3, 4, 5, 6), detto primo coefficiente di canale noto di ascissa A b il 
coefficiente C(6)=Ci(6)+Cq(6) corrispondente all'ultimo simbolo pilota (6) dello slot 
corrente (L), detto secondo coefficiente di canale noto di ascissa B e il coefficiente 
C(10)=Ci(10)+Cq(10) corrispondente al primo simbolo pilota (10) di uno slot 

15 successivo (L+l), e detto terzo coef^o' ; • canale rioto di ascissa A-l & ii* 
coefficiente C(5)=Ci(5)+Cq(5) conispt.^vH^ t fi^nultimo simbolo pilota (5) dello slot** 
corrente (L), ed il calcolo dei coefficient di canale C(k)=Q(k)+CQ(k), con k=7..9, 
avviene secondo la sequenza: 

C,(8)=[Ci(6)+C,(10^ . . . . 

20 C I (7)=[Ci(6)+C,(8)]/2 ; Cq(7)=[Cq(6)+C q (8)]/2 ; 

C,(9)=[C l (8)+C,(I0)]/2 ; C Q (9)=[C Q (8)+C Q (tO)]/2 . 

8. Metodo secondo la rivendicazione 2, in cui ciascuno slot contiene otto simboli 
pilota (0, 1, 2, 3, 4, 5, 6, 7), detto primo coefficiente di canale noto di ascissa A e il 
coefficiente C(7)=Q(7)+Cq(7) corrispondente all'ultimo simbolo pilota (7) dello slot 

25 corrente (L), detto secondo coefficiente di canale noto di ascissa B e il coefficiente 
C(10)=Ci(10)+Cq(10) corrispondente al primo simbolo pilota (10) di uno slot 
successivo (L+l), e detto terzo coefficiente di canale noto di ascissa A-l e il 
coefficiente C(6)=Ci(6)+Cq(6) corrispondente al penultimo simbolo pilota (6) dello slot 
corrente (L), ed il calcolo dei coefficient di canale C(k)=Ci(k)+C Q (k) f con k=8 t 9, 

30 avviene secondo la sequenza: 

C,(8)=[C 1 (6)+C,(10)]/2 ; Cq(8)=[C o (6)+Cq(I0)]/2 ; 
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C,(9)=rC(8)+C,(10)]/2 ; Cq(9)=[Cq(8)+C q (10)]/2 . 
9. Metodo secondo la rivendicazione 1, in cui detti coefficienti di canale da 
calcolare sono compresi fra un primo coefficiente di canale noto, di ascissa A, 
corrispondente ad un ultimo simbolo pilota di uno slot corrente (L) ed un secondo 
5 coefficiente di canale noto, di ascissa B, corrispondente ad un primo simbolo pilota di 
uno slot (L+l) successivo a detto slot corrente, essendo inoltre noto un terzo 
coefficiente di canale, di ascissa B+l , alia destra di detto primo coefficiente di canale di 
ascissa B, ed il calcolo di delti coefficienti di canale avviene mediante i seguenti passi: 

a) applicare ripetutamente in modo ricorsivo detto algoritmo di interpolazione 
10 nell'intervallo delimitate da detti coefficienti di canale noti di ascissa A e B, effettuando 

una prima iterazione in cui viene calcolato un primo coefficiente intermedio (di ascissa 
C), ed effettuando successive iterazioni dello stesso algoritmo in sotto-intervalli di volta 
in volta delimitati a destra da detto coefficiente di canale noto di ascissa B ed a sinistra 
dal coefficiente intermedio calcolato nella precedente iterazione, fino a quando non b 
1 5 stato raggiunto e calcolato il punto di ascissa B-l ; 

b) cercare, p>« ■■<■>< .••-r^ per ascisse decrescenti, un primo punto ancora da calcolare alia 
sinistra dell'iC \ v / ^fficiente intermedio calcolato, defmire come estremi di un nuovo. 
intervallo di applicazione di detto algoritmo di interpolazione rispettivamente il primo 
punto noto a sinistra ed il primo punto noto a destra rispetto a detto punto ancora da 

20 calcolare, ed applicare nuovamente, in modo ricorsivo, detto algoritmo di interpolazione 
su detto nuovo intervallo effettuando successive iterazioni dello stesso algoritmo in 
sotto-intervalli di volta in volta delimitati dall'estremo destro di detto nuovo intervallo e 
da un estremo sinistra costituito dal coefficiente intermedio calcolato nella precedente 
iterazione, fino a quando non e stato raggiunto e calcolato il punto immediatamente 

25 adiacente all'estremo destro di detto nuovo intervallo; 

c) ripetere il passo b) fino a quando non b stato calcolato il coefficiente di canale 
associato al valore di ascissa A+1. 

10. Metodo secondo la rivendicazione 1, in cui detti coefficienti di canale da 
calcolare sono compresi fra due coefficienti di canale noti a sinistra, corrispondenti agli 
30 ultimi due simboli pilota di uno slot corrente (L), e due coefficienti di canale noti a 
destra, corrispondenti ai primi due simboli pilota di uno slot (L+l) successivo a detto 
slot corrente, ed il calcolo di detti coefficienti di canale avviene applicando una prima 
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volta detto algoritmo di interpolazione per calcolare un coefficiente intermedio, 
suddividendo cosi in due sotto-intervalli l'intervallo compreso fra detti coefficienti di 
canale noti a sinistra e detti coefficienti di canale noti a destra, ed applicando 
successivamente, in parallelo su entrambi detti sotto-intervalli, detto algoritmo di 
5 interpolazione per calcolare i restanti coefficienti di canale compresi in ciascuno di detti 
sotto-intervalli. 

1 1 . Metodo secondo la rivendicazione 1 , in cui almeno un punto noto di detto primo 
o secondo estremo e un punto che & stato ottenuto mediante combinazione lineare di 
coefficienti di canale noti. 

10 12. Metodo secondo una qualsiasi delle rivendicazioni precedenti, in cui detta rete di 
comunicazione e una rete di comunicazione cellulare di tipo UMTS. 

13. Dispositivo per la stima della funzione di trasferimento di un canale di 
trasmissione in un sistema di ricezione per una reti di comunicazione, comprendente: 

- mezzi di memoria (100) atti a memorizzare coefficienti di canale corrispondenti ad 
I? uno slot corrente (L) ed almeno un coefficiente di canale corrispondente ad uno slot 

successivo (L+l) a detto slot corrente (L); 

- mezzi di interpolazione (104, 106, 108, 110) atti a leggere da detti mezzi di memoria 
(100) un primo ed un secondo operando, corrispondenti a coefficienti di canale noti, ed 
a scrivere in detti mezzi di memoria (100) un valore corrispondente alia media 

20 aritmetica fra detti primo e secondo operando, detto valore essendo corrispondente ad 
un nuovo coefficiente di canale; 

- una unila logica di controllo (102), per indirizzare in lettura e in scrittura (R/W) detti 
mezzi di memoria (100) e per controllare detti mezzi di interpolazione (104, 106) in 
modo da effettuare, mediante singole operazioni di interpolazione, il calcolo e la 

25 memorizzazione in detti mezzi di memoria (100) di singoli coefficienti di canale; 

caraiterizzato dal fatto che detta unite logica di controllo (102) esegue una serie di 
operazioni di interpolazione in accordo con il metodo descritto in una delle 
rivendicazioni da 1 a 12. 

14. Stazione radio base, del tipo comprendente un ricevitore a rastrello per la 
30 ricezione di segnali provenienti da terminali mobili, provvisto di un dispositivo per la 

stima della funzione di trasferimento di un canale di trasmissione mediante il calcolo di 
una plurality di coefficienti di canale, caratterizzata dal fatto che detta stima della 
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funzione di trasferimento viene fatta in accordo con il metodo descritto in una delle 
rivendicazioni da 1 a 12. 

1 5. Terminate mobile, del tipo comprendente un ricevitore per la ricezione di segnali 
provenienti da una stazione radio base, provvisto di un dispositivo per la stima della 
5 funzione di trasferimento di un canale di trasmissione mediante il calcolo di una 
pluralita di coefficienti di canale, caratterizzato dal fatto che detta slima della funzione 
di trasferimento viene fatta in accordo con il metodo descritto in una delle 
rivendicazioni da 1 a 12. 
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TTTOLO 

Metodo e dispositivo per la stima della funzione di trasferimento di un canale di 
trasmissione 

RIASSUNTO 

5 Un metodo per la stima della funzione di trasferimento di un canale di trasmissione in 
un sistema di ricezione di tipo UMTS prevede il calcolo di una plurality di coefficienti 
di canale, compresi fra coefficienti di canale noti corrispondenti a simboli pilota, 
mediante la reiterazione di un algoritmo di interpolazione atto a calcolare un punto 
intermedio (Z, f(Z)) fra un primo estremo ed un secondo estremo di un determinato 

10 intervallo, il primo estremo essendo costituito da almeno due punti noti ed il secondo 
estremo essendo costituito da almeno un punto noto, il punto intermedio che si vuole 
calcolare avendo per ascissa (Z) il valore di ascissa del punto medio fra i punti che 
delimitano Tintervallo arrotondato all'intero pid vicino al primo estremo, ed avendo per 
ordinata (f(Z)) la media aritmetica fra Tordinata del punto noto del secondo estremo e 

15 I'ordinata di un punto, scelto fra i due punti r • - v Imo estremo, avente una distanza* 
dal punto intermedio pari alia distanza f*a\Tt ; feitermedio ed il punto noto del^ 
secondo estremo. 
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